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ams thesis  anelw@es= 2 survey of ane literature througn 
1967 on the methods of determination of the thermodynamic proper- 
ties of nitrogen-oxygen mixtures in the two-phase, liquid-vapor 
ae 2 OY) “Tne virial expansion, empirical equations of state, pseu- 
mocritical Bere rants excess functions, latent neats, and grapni- 
ewemetvnods are considered, 

The methods just enumerated are evaluated for their simpli- 
Canty , eee and accuracy. On the basis of Une stale sone 
ere and the: availability Of a2pDrerriave wate, thermodynamic dia- 
Meemceeare conStructed for 1, 5, 10 and 20 atmospheres of pressure 
in the two-phase region. Diagram coordinates are enthalpy- 
composition and entropy-composition. The diagrams at atmospheric 
pressure show good correlation with known data, ‘while those at 
higher pressures are best estimates using excess functions and 
ier! iZing assumptions. 
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LIST OF SYMBOLS 


+ x oh % 
Oe Ons 49 bj;~~~- Single substance, mixture and pure com- 
ponent of mixture constants in Hedlich- 
KWong €On. of —Staven 


a, A, b, B, c, C --- constants in Barker's expression for Gc". 


Eur), C(T) ~-~ 2nd & 3rd virial coefficients for a sinmeleucin— 
stance, 


Baye Ee ex) Ba , 
~-~- 2nd & 3rd virial coefficients for a mixture. 


Cy ake C (ieee. Ci. 


Bane Bape Crane CaRR --- 2nd & 3rd virial coefficients for pure 
components in a binary mixture, 


Bape Capp CaapB --- ena & 3rd interaction or cross Virie lecoeci ae 


clients for & binary mixture, 
x x 
mie, C(T) -~- 2nd & 3rd universal, single substance, virial 
COcil Jouenvs. 
b, --- "hard sphere" molecular volume =(2/3) NO ? | 


Ce -~- Specific heat at constant pressure for perfect gas or at 
zero pressure, 


Mig) etc, —--- 4th virial coefficient. 

E --- internal energy. 

f --- fugacity of a single substance, 

i --- fugacity of component i in solution 

oe W-~= fugacity of pure component i at same T,P of the mixture, 
G --- Gibbs free energy = H - TS, 

H =<--~ enthalpy = BE + PV 


a [AB ih differential heat of yaporization and condensa- 
Y tion at constant temperature. 


j --- absolute joules 
k ~~~ Boltzmann constant = 1.38054 x 10°? j/oK 
L, M, N #-- defined by Eq. (58). 


n --—- number:of moles. 


LV 





pene 


Peoolit. 


Ps ---~ pure-component pressure of component i in @ gas mixture, 


P, -- 
>) 
) —_ 
a 
: = 


ee 


Xyoo- 


ee i 


Super 


ae a 
Les 
a 
i = << 


O eo aa =p 


- partial pressure of component i in ideal gas mixture, 
Pe = xy PM Dakton’s: lanje 
universal gas constant = 0,0820535 liter atm/mole °K 
any extensive property 
Dawidaloapropertby of component, 1 in solution, 
entropy 
absolute temperature, 
~- base or reference temperature, 
intermolecular interaction energy. 
molal or specifis volume 
- partial molal volume of component i in solution, 
mole fraction of a component i in a mixture = n,/n. 
compressibility factor = PV/RT 
‘Soneayenas 
EeecCos Properly. 
property of the liquid phase; property of the liquid line. 
property of the vapor phase; property of the vapor line. 
BPeQvoertLy Of mixing. 


zero pressure state or perfect gas state at any pressure, 


subscripts 


b --- 


1 e es aa 


m en Ge G= 


CG oe a= az 


on 


refers to the normal boiling point. 
pure. Component property of a mixvure, 
Paoperuy of & mixture. 

Ceibacal property. 


- pseudocritical property. 





Yr e=-= reduced state, 
S -8-- reference state. 


Greek symbols 





X --- volatility ratio, 

 —— AGrmivily COcinicicny, 

a -~-- general function of temperature. 
ee--—- characteristic molecular energy. 
a characteristic melecular distance, 
es frm ve change Of Quanvircy. 

[U--- chemical potential = G 


mmeescacte properties are given in their molal or Specific proper-—= 


ty form unless otherwiss noted, 
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It is the purpose of this thesis to present the results of a 
search of the literature and the Subsequent analysis of the find- 
anes of that search. The subject of the inauiry was (Uneevacnne] 
aynamic properties of nitrogsen-oxygen mixtures in the two-phase, 
liquid-vapor region. The reasons for the study of the elemental 
properties of any basic substances are self-evident especially 
her substances as fundamentally important as the mixtures of nitro= 
sen and oxygen. It should suffice to say that the nunt carried 
the author through the fields of cryogenics, low temperature phy- 
fees, Nigh pressure physics, thermodynamics, physical chemistry, 
Memecvlar chemistry and statistical mechanics, It 1S the intent 
fois LUnesis to extract from these areas of science and engi- 
neering the most appropriate analytical, graphical and machine 
methods @f determination now available and to present tha pre 2 = 
ieeoeaetermined therefrom, 

Acknowledgement and appreciation are expressed for the 

Pea, assistance, and direction of Professor Ml Pete Jee 


meer ne MIP Cryogenics Laboratory, 





The direction of this work consisted of two distincive a2. 
The first part of the effort was directed toward a thorough sSearcn 
of the scientific and enginecring literature. This Wastecce — 
plished by a personal search of the Chemical Abstracts through 
1967. Machine programmed searches were conducted for the author 
by the Cryogenics Engineering Laboratory of the National Bureau of 
pramaards and the Defense Documentation Center. These initial 
Gii1oOrtcs were principally centered on nitrogen-oxygen mixtures and 
the thermodynamic approaches to the derivation of properties of 
binary mixtures, As the work progressed, molecular theory became 
of interest, and the extensive bibliographies in the major texts 
Saeeme ti1e€id were utilized. At the end of each subsection of the 
Results, the references are listed which deal explicitly with 
tMecu Subsection as annotated in the text, 

The second part of the program was the analysis of the data 
Which had been found. It must be said that the determination of 
the properties of a mixture would be a relatively easy problem if 
mieemixture behaved Simply as the average of its constituents, 
Maportunately, not only are the constituents themselves not ideal, 
bueeetiso they interact with each other in a non-ideal manner, 

It is intended in this thesis to discuss the PVT behavior of 
mixtures and the various methods of deriving théir tnermodynamic 
properties: enthalpy and entropy. Specifically the next section 
will cover the following: 

a, the theory of corresponding states and pseudocritical con- 


Steancts. 





b. empirical equations of state, 
e. the virial expansion, 
Gee partial property correlations and mole fraction averaging, 
e., excess functions as deviations from ideality. 
Moe eraphical derivations. 
eam Jatent heat techniques, 
No attempt has been made to make interpretations based on statis- 
tical mechanics and its application to molecular theory, The read- 
Pmerncverested in approaches based on this field must consult the 
works of Eyring, Kirkwood, Lewinson, Hirschfelder et al for their 
Greatment on the formal and model approaches in this area, 

The final section of the Results will correlate the discussed 
data and provide a Saar! presentation of the enthalpy and en- 
tropy of binary mixtures of nitrogen and oxygen at various pres- 


Sures, 





tf 


HeoU le ie 
Phase Fauilibria 

panies the information most consistently reauired in calcu- 
lations concerning the two-phase region is the equilibrium behav- 
jor of the system as defined by Gibbs [2]. The study of vapor- 
liquid equilibrium would obviously be simplified if the phases 
behaved as ideal solutions, for then we should only require infor- 
Beeson about the pure components in order to develop yalues for 
ie mixtures, 

Deviations from ideal behavior occur in both the diauldgadd 
the vapor phases as & result of the forces that exist between the 
molecules, The magnitude of the intermolecular forces depends on 
the nature of the molecules and the average distance between them 


as expressed by 


Usn(t)/E ap = f(0/C ap) (1) 


This expression of the principle of corresponding states [2,3] 
says that the interaction energy of a pair of molecules A and B 
at a distance r apart may be expressed in terms of a universal 
moinewton £ for all A and B, The parameters C ap and O an are a 
characteristic energy and distance respectively for the pair of 
molecules, In the liquid state the molecules are relatively 
close together, and the force fields are strong, so that differ- 
ences in shave and charge distribution will cause deviations in 
behavior, In the gas phase the force fields are not as strong, 
and at low pressures ideal mixture behavior is a good approxima- 


tion. Note that a perfect gas is defined in terms of negligible 





forces between molecules, whereas an ideal mixture or solution 
requires oniy the same forces between like and unlike molecules. 
Hence, at low pressures ideal solution deviations are Drine? pale 
due to the liquid phase while at high pressures both vapor and 
liquid phases contribute, 

It is perhaps instructive to review the phase behavior of 
mixtures under pressure [4,546]. In Fig, 1 the curves OCy»o and 
OC oo represent the usual vapor pressure curves for the pure com- 
ponents terminating at their respective critical Points. | For wome 


binary system, however, the phase boundary is the region enclosed 
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Figure 1. Phase boundaries of pure components and mixtures. 
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ie ene  Cuky ec OC_O within which both the liquid and the vapor COeCX~ 
iet.  Fhescritical point.of the mixtimme. oa is defined as the 
state at wnich the properties of the vapor and liquid prasesoma 
equilibrium are identical. If the pnase boundary curves of Severe 
mem xtures of different compositions are drawn, te vocircecr 
fe@eir critical points form an envelope indicated by Une vg2sed 
line in Fig. 1. The maximum pressure and temperature for the ex- 


merence of twoephases= 2re~ cakled thes crmicondembar, C and the cri- 


Pp? 


Gondenthernm, Cis BeSPeCC UI Ve Ly 17]. These are illustrated in Fig, 


2as is the phenomena of retrograde condensation 8 , 


Cri cacal 
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ioe, 2 it can be seen that if a line of constant UCémperarure 
meetollowed, condensation will begin at point a and will continue 
eh decreasing pressure until maximum condensation is reached at 
point b, at which time vaporization will occur until the system is 
again all vapor at point c, A similar phemomena will take place 


meeconstant pressure witn vaporization occurring with decreasing 





vemperabure Trom pDoOINUvS a Uo e. 

The diagrams most frequently used in representing ecu ieae 
um data are those of pressure or temperature as a function of 
molar composition. Fig. 3 shows the trend of the phase boundary 
curve with increasing pressure when the mixture critical tempera- 
ture and pressure fall between those of the pure components, The 
area of two-phase coexistence decreases with increasing pressure, 
The pnase boundary curve breaks away from the pure component whose 
Critical temperature is lowest and recedes along the critical en- 
velope, finally disappearing altogether at the higher critical 
pPoinc. 

Turning to nitrogen-oxygen mixtures in particular, the earli- 
Pomeinvestigzation of their equilibria was’ conducted by Baly in 
1900 at atmospheric pressure [9]. The temperature scale used by 
him was in error, and his vapor composition measurements proved 
wembpe innaccurate, Inglis /10] in 1906 made measurements at 
2b .9 and 79.07°K resulting in vapor pressures between 1 and 32 
atmospneres, From this work an empirical formulstion was made 
for the tabulation of isotherms and isobars,. 

During the next few decades the only measurements were those 
of Trapeznikova and Shubnikov [13] in 1934 on an isotherm of 
85°K, and the static measurenents of Sagenkahn and Fink [ 14 at 
l and 2 atm in 1944, 

iitcle else was done until 1955 when Armstrong 2nd aise cen- 
Meaenes at the National Bureau of Standards used a circulation 
method to make accurate determinations of the compositions and 
vapor pressures at 65, 70, and 7 ao [15]. Two years later 


Cockett conducted extensive measurenents at 1.3158 atmosovheres 





[16] with ‘results reasonably close to those of DodSse and vunee 

DeCent. <i 2a 1eant fess Senn are Sleime Gtt Din [17 | who 
Peed the sraticemetrod for nine intermediate mixtures at pressupes 
to 10 atm. His results likewise show those of Dodge and Dunbar 
to be reasonable as do the experiments of Latimer [18 | in loo7e 
In 1961 Yendall and Olzewski followed with a report on the satur- 
ation properties of nitrogen-oxygen mixtures hi} . 

The WADD Technical Report of s1o6u [20 | presents a compila- 
tion of the data of Armstrong et al, Dodge and Dunbar and Sagsen- 
kahn and Fink on T-x eee The data of Sagenkahn and Fink at 
ieana 2 2tm and of Cockett at 1.3158 atm are presented on P-x 
diagrams. The graphics use a base of Hilsenrath's | 21] pure com- 
ponent vapor pressure values. Tables of selected experimental 
values of all of the above are given. | 

hir Product and Chemicals, Inc. also publisned wndér govern— 
Mm oeeconvract a technical report | 22 | with a P-T diagram of vapor 
‘and liquid curves based on the data of Dodge and Dunbar for the 
two-phase region and the vapor pressure data of Hoge for oxygen 
[23 and ‘Michels ey ale hrevedman saad white, Dodge and Davis, and 
Henning for nitrogen [ 24,25,26,27]. The critical envelope is 
based on the work of Kuenen et al [ 28,29] and can be considered 
peceise to LOE. The work of Kuenen et al also indicates the 
presence of retrograde condensation over small ranges Near cine 
Smeeercal voints. , 

The most recent work is that of Haselden and Thorogood in 
1963 at high oxygen concentrations | 30 | which showed a good corre- 
lation with the earlier works of Din and Latimer below 96 Zmole 


oxygen. They demonstrated, however, that the relative Vola U7] Tae, 


Malls off at high oxygen concentrations. This is alllered to te 





indicative of the occurrence of association, that is, the weak 
chemical interaction between the oxygen and nitrogen molecimecs 

For tne work at hand, reliance has been placed primarily in 
the data of Dodge and Dunbar, because it is considered to be con- 
feeeucnt and detailed over a wide range of values ee bieciee ian. 
can be considered accurate to 10.1°K at pressures near atmospheric 
and to about 7% at pressures above seven atmospneres [22]. 
saturation values for the pure components have been taken for 
This work from Stewart for oxygen | 31] and Strobridge for mavre— 
gen [32] , Since the data of these investigators was used for all 
other pure component property values. Stewart's empirical vapor 
Pressure equation was fitted to the data of Hoge | 23] while that 
of Strobridge fits the data of Friedman and White | 25| and Arm-~ 
strong [33]. 

Modividual temperature-composition diagrams heave been drawn 


moet, 5, 10 and 20 atm and are presented in Appendix A. 
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ioe. Vitel 2eue clon 
The virial equation of state is simply 2 Series expane io men 


me compressibility factor in termssonudere 


Z=PV=1, BIT x) , C(T»x) , DIT x) | | (2) 


RT V (emis OCC , 


mecause of the molecular interactions, the virial cocifticienrees2ne 


functions of both temperature and composition and can be written 
as [a | 


B B(T,x) , Bes x and (3) 
= “ED 


©2 
lf 


2 n 
. c(T,x) =). 2, 
‘i=A j=A 


For a binary mixture these become 


ire 


Xs x 52K C4 jk | (4) 


2 2 
B= XaBaa + .2XaXBBap + XpBpp and | (58) 


C 


i 


: 2 2 3 
m= *acaka + 2XaX%pCaap + 3X,%pCppa + XBCoRap (6) 


where X, 18S the mole fraction of component Pes 

Molecular theory tells us that for non-polar molecules (angle 
independent molecules possessing no permanent charge distribution) 
such as those being investigated, calcusation of their virial co- 
ere rcients devends on the chojce of proper potential functions. 
The potential function represents the attractive and repulsive 
forces between molecules and many are found in the literature 


| 2, 3]. Although these forces are generally exponential in char- 





acter, they have generally been simplified to reciprocal power 


morms « 


The most popular relationship is of the form 


U(r) = K rvs Kor7™ Ce) 


where the first term represents the repulsive potential energy 
ome che second that of the attractive contribution. Where yn and 
m are chosen as 12 and 6 respectively there results the well known 


Lennard-Jones potential (compare with Eq, 1). 
U(r) = 4E (G/r)1*-( G/r)® (8) 


It can be shown 4 | that by using the Lennard-Jones potential the 
meen component virial coefficients can be expressed as 


a 
S 


B(T) = f(T 4b.) or B (PT) = B(T)/by= f(T") (9) 


where be is the "hard sphere" volume defined by (2/3) 7C! 103 and 


3- 


omry(C /k). Similarly 
Bey = C(T)/b,? = £(T*/€ /x) 7 (10) 


Bas. (9) & (10) nave both been tabulated and plotted. 
ge interaction or cross virial coefficients are .sicae a 


woe THOre Gifficult to determine. For an ideal solution obey-= 


2 
Caap=(2CanatCanp)/3 - (Baa-Bpp) | (11) 


2 
Capp=(Cannt2Copp)/3 - (Baa-Bpg) 





mogen its the Lewis-Handall rule for mixtures of simple g2seq and 


mesilts in the followings for a binayyw as tire. 


By =X, Ba xtXy8ap and 


2 
Cm=xaCnantXpCpyp - 3X,Xp(Baa-Bap) (12) 


? 
It has been shown [6], however, that as NN COU ae Le @ 


universal function of Stl Swe $0. G00 is Bp oe) 1G) yee havens wall 
ony KT/€ ape The measurement of second Virials i uem cs Oneal 
lute gases where the third virials are negligible can be used to 

Mereiy the various combining rules derived on the basis of This 


relation. The most popular and simple of these are the Lorentz- 


Berthelot rules | 7 | 


Eape(€ ar€ pp)? and Oap=(O ant Opp) /2 (13) 


which is based on London's theory of dispersion energy [8,9,10], 
Meese rules give a better account of the second virial coeffi- 
cient of simple gases than the Lewis-Randall rule, yet there is 
Meal doubt || Waetner Any rule will yield accurate valneewen Bap 
emer) under tne best conditions of Seale non-polar and spherically 
Symmetric molecules, To make matters worse, very few experimen- 
tal measurements of second virial coefficients have been made on 
mixtures of simple molecules. Newer and more accurate functions 
such as the Kihera povential are presently coming@wyncer more aud 
more study [11,22,13,14]. 

men CONSIGeration 15 On Mixtures otner erameq tile w2esces, 
memmexect relation is available to correlate the interaction ener- 


fies, Present practice | 6 | 1S to extend thewemame ole, by 





considering an average interaction energy 
pee )> / © Gans | r/ O(x)| a) 


Magie the virial equation provides the best theoretical basis for 
Beapimating the PVT properties of mixtures, 15 cannot deal satis- 
meacecvOrily with any but Low amecs oes FAS Mix GUC Seem ae eee eae 
femmes uimate third and nigher virial coefficients genecrarves an 
Overall loss of accuracy, 

The thermodynamic proverties of enthalpy and entropy are de- 
mavable from the virial equation [ +,15,16,17 | . PROM Ehe teenar 
and Second Laws for a single-component, closed system dE=TdS - PdVv 
mien —s+PV whence dH=TdS+VdP, If this eee 1S divided Dyan ard 
fis held constant, the relation (OH/OP)p,=V-T( OV/ OT)p results. 


If Z=PV/RT, then ( OH/ OP) p=(-RT“/P) (0 2/0 T)p or after integration 


P2 
(Ha-Ha Iga-R04 (32/0 T)pd(in P) (15) 


fer, —o end T,=T/T,, Pr=P/P,, then 


Py 
(HAH) /atgann? | ( d2/ dT ,)p d(InPp) (16) 


where (H-H?) is the enthalpy change from P=0 to P=P at constant 


termperature. Then for mixtures Git Eas. (2) & (15 )eawe 
0 @ ee 8 
(H,-#°) /at=1/v [ 3,-ras,/at| + 1/v?{ c,-T/2(ac,/aT) | + (17) 


From the Maxwell relation CO.S/8O) = (OV/OT)p and Z=PV/RT we 
have (OS/O0 P)p=-R/P [ Z+T(D2/DT)p | . Since the entropy blows 


Mees pressure goes to zero, it 1S convenient to use the artifice 


Peeincvezrating this relation as @ real gas to a pressure near 





Zero Where propertics approximate tneose of 4 perieccne eae 
equation is then integrated using Dewpecc cas proper a. tou 


eonvenient reference pressure, PF Se thay Low pressure wmc cre 


q! 


on the entropy mutually cancel Be Then 
Ie ie 

(s-S,)/ref (1-Z)d(1nP)+1in(P_/P) = af (OZ/od7T). d(inP) C16) 
0 : 0 P 


or in reduced terms 


Je 7 P 
x x : 
(8-8) /A-f (GFZ) aeNa) +n ae a vf (02/ OT, )p a(inP,,) 
(19) 
where (S-S.)/R is the entropy change from 


ie. to T,P or T.,P,., to T,,P,. Then, Eq. (2) we Nel, Vids) feaiice 


2 2 2 
(Sp-S,,¢)/R=-1nP,~T/v| ap, /at | - Be /2v"4,,/2V“ - r/2vic,, /ar | (20) 


Values can be determined for dB/dT by differentiating curves of 
B¥(T) as a function of T, using the relation dB /aT =(€ /k/b,) 
aB/dtT. 


mone liy to determine An and Oe theo re lagioenc 


O =O" : | | 
Haz Lx hy and Se 2 X5556 - R Linx, + RinPs (2) 
Gan be used, where Hy is the enthalpy of component i at P=0 and 


oa Mmenene entropy of component i at reference pressure, Poe 
The virial equation is not the most convenient method of deter- 
mins the thermodynamic properties of ‘mixtures unless Comet: 
facilities are available. If a virial equation has been formu- 


lated for PVT properties, however, it would not require much 


metry ionel effort to program the derivatives of the equation, 
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Pseudocritical Properties 


It was stated earlier that the molecular energies of inter- 


action could be described by their dependencies on @ universal 





ay 


function as expressed by Eq. (14). The molecular theories of 
Prigcosine et al [3 | in which explicitGepele ven ce tors Une se mene 
ters € (x) and CG (x) are developed in terms of their pure and 

ant © nn» Gp, rr through the 
Lennard-Jones or other m-n potential forms, are formally equival- 


mixture forms, ie, Cie € 


eae to the pseudocritical methods of Kay eft al Ee Mie  selataens 
KT/ € (x) and V/G 7(x) are equivalent to universal functions nae 
may be replaced by T/T (x) and V/V (x), where T Ax) and V (x) 

@re pseudocritical properties of the mixture [3]. The pseudo- 
Beau1cal concept is based on the existence of a mythical pure 
Mmueicd 2t the termperature and pressure of the mixture which exhi- 
bits the same properti?s as the mixture, The critical properties 
Ore cnis fluid are the pseudocritical properties of the mixtures, 

A sizable amount of literature has been written on the deter- 
mination and use of the pseudocritical constants [ 241, 5,6,7]. 
Kay's simple mole fraction averaging rules | 2| are The mOsiteatnre. 
Milently used; in fact, the various rules for pseudocritical temp- 
Seepure all reduce under various Simplifications to that of Kay, 


ley pee 1 \X)> Kay's rules are given by 


7 2 we (22) and P ).x,P., (23) 


Cm i ab Cm i Cl 


For T,,: Kay's rule produces values within 2% of other rules for 
05s To /T op Z20 and 0.5= Pa,/PoR <2 [8]. For Valuesmo. P adif- 
fering by 20% or more, the modified Prausnitz and Gunn rule | 9 | 


Should be used, 


a i eee 
em ee ven en = ; ea ; xT 47 7 Oa 5 (24) 





ray Ss rule has been verified to Dredic gar nc true CriGiee lav emi. 
erature of NA7P9 Pee ies. bo witnin Oia 

Using the pseudocritical rules (opie tres Olen b One ead 
mayeen, the author computed enthalpy deviations by the method auce- 
Rested by Heid & Sherwood. A Sune Of “Girss Me tiicd= follows: 


1. Pseudocritical constants for the mixture were calculated from 


the following pDure component values: 


Nitrogen Oxy cen 


iam 126,26°K 11,14] 154.77°K [12,13] 
P. 33.54 atm [11,14 50,14 atm [12,13] 
v. 90.07 em?/mole [15] ,075 em?/mole [ 12] 
Z 291 [8 | .292 [8| 


2, Base enthalpies, at 


based on the enthalpies of the perfect gas at 1 atm at the 


for nitrogen and oxygen were selected 


Memmal bOoOlling point and triple point respectively [11,22]. 
3. The ideal state enthalpies at the various isotherms of inter- 
est, see are then calculated from 


? | 
O O O 
Hn, = Hoff C., aT (25) 

O 


where CoNeT 29 +121 j/mole°K and C oa 29.300 j/mole°K [16]. 


4, At a perticular tempereture the perfect gas mixture enthal- 


p 


se O 
pies, HH, may be calculated from 
m 


4 


oO. ) _,0 
mee / xi, (26) 


i 


5. The mixture enthalpy at the system pressure is then 





eee epaed)) hae te | ee.) 


where the isothermal enthalvy deviation En ~ Ho can be deter- 

mined from Suitable tables [17] . 
Reid and Sherwood predict errors of 10% in the vapor phase and 
20% in the liquid for this method [8 |. This author carried out 
calculations at 84°K and 1 atm and at 112°K and 10 atm. In the. 
Pepor phase for the case of 84°K and l atm, although the calcula- 
ted values differed from those obtained using a mole fraction 
average of the pure component enthalpies by as much as 35 j/mole, 
the results were linear and had the same slope (Fig, 4a). In the 
mauid region at this temperature the results showed a small neg- 
evive deviation from linearity (experimental evidence indicates a 
positive excess heat) with the absolute difference increasing 
with increasing %mole Oo (Vie. tb), At the nisher recaineuine 
calculated vapor isotherm was no longer linear, but exhibited a 
pazecable negative excess heat. Likewise the liquid isotherm at 
this pressure displays a somewhat larger contrary behavior than 
before. All isotherms are extrapolated into the wet vapor region, 
These calculations seemed to be Fadlesiiemie for their Geviarious 
from known behavior, ie, negative excess heat. The tables offer, 
however, a relatively simple computational technique and are ap- 
plicable to the liquid phase as well as the vapor, Perhaps better 
values of the pSeudocritical properties would produce more logi- 
Cal results, but for the present, Simple mole fraction averaging 
Seems to provide more satisfactory answers, 

While Zabiles Sieh as tUnose Of LyGderson €t-ale doenee to ao wore 


entropy deviations from the perfect gas base, a relationship can 
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Figure 4a. 
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Figure 4b. 
Figure 4a & 4b. Liquid end vapor lines determined by various 


methods. Dashed lines join pure component values. 





be written between properties of other values which are tabulated, 


aa (18) can be written @ 


ff - S.)/R=(H - Hoy Rin shia Gay eres (28) 


when P=. and yates "of both (i =) Ho ecm Ci ioe) have been tabu- 
faved. 
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Empirical tauations of State 


Perhaps the most underdeveloped of the techniques for the 
mevermination of the properties of mixtures iS that of a workable 
eemavion of state, particularly in the two-phase region. 

The Benedict-“Webb-Rubin equation [| witn eight constants 
meemoeen fruitfully tested, but unfortunately is @pplicable only 
wemiebeznic hydrocarbon mixtures, 

Some effort has been made to determine coefficients for the 


Peo constant Redlich-Kwong equation [2] as follows: 
» + a/(T°*Fv(v+8 ]| (V - >) = RT (29) 


+ % 4 % 
Giieme fOr &2 mixture a = .: AM Oe = eS 
mM 7 ah) 2a m 5 se | 


x? 2.5 
no 
ard also a, = see Pai 


an 2 


It is revorted that some success has been obtained in using 
meauwced temperatures and pressures derived from Kay's rules in 


the generalized Beattie-Bridgman equation [3]. 


i ea 


ee (T /V 5) (2 : 0.05/TAV 4) y a 01867 (2 : 0,03833)/7,4]} 


-(0.4758/V,.2) [ (2 = 0.1127)/V,,, | (30) 


Recent efforts [4,5] nave been fairly Successful in treating 


Bas mixtures of CO. and No With mixture rules fortene ter v1 ou 


2 
Sewation of state, 


5 
Pe 3 f.(vV -»b)7* where f, = A, + B,T + C,exp(k'T/T,) (31) 
mee, Ct : - ° 
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None of the equations of state have been extensively tested with 
mexcures, With tne posSible excenvtion On seo ioc ola 1 oo 
ieeent hydrocarbons, and all require machines nee eae 
derivation of thermodynamic properties using #as. (15) & (18) 
the results will be only as good as the equation used, and some- 


memes worse, because of the errors involved in the determination 


of the derivatives of the equation, 


Bererences for the Empirical Equations of State 
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Partial Provertvies 

It was mentioned earlier that use mignt be made in some cases 
er partial properties in determining mixture properties. Partial 
Meee Droperties define the properties of components in @ mixture 
Pewooposed to those of the pure components and represent a part 
Memes LOLA) property of the particular mixture. Im particular, 


es property is defined as 


Ry = WOR/ OTe, Ps or K = di nk (2) 


where R is any extensive property and n, is the number of moles 


eeecomponent i in tne mixture. These relations are governed by 





mae Gibbs=-Dunem equation 


yy ny GR; =- Q esis) 
i 
men = R/n = e x,H,, then it can be shown /2 | that 
i 
oe & - x, ( OR/ See) eile = OO) | (34) 


The Lewis relation [2 | for fugacity agit Ys prelateduroutac 


chemical potential can be expressed for a mixture as 


M, = (OG/ Oxy PyT = RT inf, + Q,(T) (55) 


where 0. (2) is some function of temperature only. As the pres- 
sure approaches zero, f approaches P. as opposed to the chemical 
Potential which blows up as can be seen from (OM/OP)4 = Vas 


At constant temperature and composition Ea, (35) can be writvcen 
[ 3.1] 


a E one 
RT In(T,/f,x;) =f ai 
5 ae . 


, ~ V,)aP SRT in, (36) 
where un We “derined 2s the activity coefficient. fer mesmo 


utions this reduces to’ the Lewis & Randall rule [2 | 
if, = ET, such that t, = 1. (ey) 


Mate in an ideal solution it is only necessary that the inter- 
merecular forces be the same for all components, in @ perfect gas 
mae molecular forces must be zero. In this case the pure compon- 


eime tugacities are equal to the vressure and Raoult's law results: 





> ae . 7 (38) 


It should be mentioned that the Lewis-Randall rule generally 
provides a good approximation for most real non-polar gas mix- 
fres at any pressure as long as the component in question is Sin 
jlar to the other component or present in excess of the other com- 
ponent 5 |. 

If the fugacity of the pure components in the mixture can be 
Secermined, the enthalpy deviation representing the non-idealbity 


of the mixture can be determined from the relation 
= | 7 O Z 
- Ld fet x,)] sar |p, = (H, - Es : (39) 


The mixture enthalpy then, follows from the partial property and 


Periect gas base state relation 


, (lee me H? | (40) 
m A 1 i i ak 
where the pure comvonent perfect gas values are determined as be- 


fore, 
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Excess Functions 


Before discussing deviations from ideality, it would do well 
morcefine somewhat more fully the ideal mixture. One definition 
[2 | states that it is a mixture whose components have chemical 


potentials given by 
ML, (P,T,x) SOE) a tnx, (41) 


where {I ,(P,T) isthe pure component povential av them pre ccm. 
and temperature of the mixture, It follows that for an ideal mix- 


ture 


>) 
if 


dx My ed dx, Inx, 
(42) 


XH 


moo, + R - Loa. 
1 i i 1 1 


co a 
ees 


‘ps [7] vol a +[~] 


iene Log sums are referred to as the ideal free energy Of mixing, 
G", and the ideal entropy of mixing (Gibbs paradox), S™, The 
ideal free energy of mixing is always neces and that of entro- 
py always positive. The ideal enthalpy and molar volume of mix- 
ing are equal to zero, 

Excess proverties may be defined in these terms as the dif- 
Mmerences between the values of the proverties in the real mixture 
and the values of the properties given by Eqs. (42) for an ideal 


mixture of the same composition. 
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It has been Seen that the chemical potential may be related 
me the activity coefficient, so that fe=ea wee 4 en are ers aol, 


pressures 


Pe (P,T,x) -{1,(P,T) = RT inx; f, (43) 
Tne partial molar properties may then be given by 


H, = 4H, - RT“(O1nf,/dT)p 


S 


ee, - = inx,/, - #10 Wx,/d7), (Uy) 


V 


p= Yt RT(D1n J 4/ OP) 


The excess thermodynamic properties may then be obtained by Sub-~ 


pracvting° the ideal quantities from the vartial molar. quantities. 


E 
G Roe ) x in 
7 5 if 


(45) 
E Z 
; | . 
pe = -R dkny 7) Se stik Ys x(01nT,/ 07); 
B 
7 = RT dix (inTy/ OP). 
maceexcess functions are related by 
Pe . tse . Ae 4 5p dc"/ dT), ; (46) 


Considerable effort has been made in the past decade in the mea- 


mmement of excess functions, In particular, regarding No-O5 


E 


TG e 
mixtures, Knobler et al [2.3.0 | measured H and V’ at Tees Wil- 


:y 
fein & Schneider [5], co” at 77.59%: Pool et al GE GY at 83.8%; 
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and Duncan &« Staveley [7], — at SG ee in addi Gicon,sveives 
of oe have been calculated by some of these authors from the 
vapor pressure measurements of Dodge & Dunbar ls], Armstrons et 
al 9}, and Din [10]. These results are plotted in Fig.) JZ or 
equimolar liquid mixtures, | 

It can be shown [2 | that the variation of vap0r pressure 
emiucans Lemperature can be related to the excess Gibbs irvece en- 
ergy. Vapor pressure and composition can be measured in an 
equilibrium still, a large number of which have been devised 
[a1-2u] . Tne principal problem in aid stills Veetne doe ae wine, 
in eliminating vapor condensation before it reaches the vapor 
moons DOint and of preventing liquid droplets from being car- 
ried over to the same point. These problems can be eliminated by 
Meenods which utilize only the total sata vapor pressure as a 
Mmemevion of liquid composition. Then the vapor comvosition may 
be found analytically using the Gibbs-~Duhem equation. Several 
methods have been used to measure the total vapor pressure | 25~ 
29]. 

If isothermal vapor pressure measurements are known, values 
of ch may be derived from knowledge of P, ee x” and the second 
virial coefficient. If the easier experimental procedure has been 
followed, which measures xt and not x”, there are several graphic, 
algebraic, and machine methods which may be used to solve for e 
|27,30-35]. These Solutions are more complicated than those used 
when both compositions have been measured. Other methods have 
Meen devised for the calculation of the liquid composition, if 
only the vapor composition has been measured | 36-38]. Finally, et 


eee Our properties are known, then the thermodynamic consistency 
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of the results may be tested with the Duhem=Margules equation [39] . 
The most widely used test for consistency is that developed 
by Coulson & Herington [uo] and by Redlich & Kister [wa]. If the 
moleatility ratio is defined by X= xp/xp / 5h where B is the 
more volatile component, then the difference in excess chemical 


botentials of two components can be given AS 
3 D it B,.2 
Ma - Me - ap int G/ Td” = (06"/Qxt), = RT Inox = EES) 


Z 
(Vx - By) - (P = Py)(Vg - Baa) 4 2P(OB)anixq = x8) (22) 


ish 
By integrating this expression over — Where G is zero at x=0 
and x=l, we obtain 


| \ 
frnco P,/P, )dxy _ var f [(P = aCe Zs Bap) Sy P.)(Va ie Baa) 
© oO 


~2P(O B), ara xv) | axt, (48) 
If the results are consistent, the volatility ratio will equal 
the pure Ole erate pressure ratio, ie, & =P L/P ys at the composi 
tion where oe is &@ maximum or @ minimum value ina non-ideal mix-= 
ture in equilibrium with a perfect gas, Therefore, by DIOt time 
KP,/P., 25.2 [unctlion of x the areas above and below the axis 
Will cancel if the results are consistent, 
Unfortunately, it has been demonstrated [28 | that tnismtess 
fails for liquid mixtures in the vicinity of their norral boiling 
points, because it assumes tnat the vapor benaves as a perfect 
&2S. Rowlinson [2 | modified the method by noting that (OB) 


AB 
made no significant contribution and that P could be peniaced by 


Raoult's law pressure, xlp + xtp SO Ena 


A A Bae 





l 
[ance P,/P, axe = iA ~ P53) (VA - Vi =a 2 Bap) | / 2a (49) 
an explicit relation which is convenient to use, 

The most widely noted method of determining G® from measured 
isothermal values of total pressure and Ii@uid (conposi tomes that 
or Barker [ vi]. This autnor has made calculations of eh for 
N5-O5 mixtures at intervals of 10°K fron 80°K to eo using 
this method. The procedure is outlined for interest in Appendis eee 

Rowlinson / | has indicated the difficulty in defining the 
excess free energy at vapor pressures greater tinea 2 to 3 atm: 
because of compressibility effects. This can be Seen imate 
where the values of cE calculated by the author have been plotted 
With those calculated by other investigators at lower tempera- 
tures. It should be noted, however, that the slope of the values 
between 100 and 120°K (dashed line) is identical to that at lower 
temperatures (solid line), and one wonders if pernaps a method 
mor calculating ce including the third virial coefficient might 
not displace these values into line, 

At any rate, studies in the vicinity of the two-phase Per ron 
at pressures higher than atmospheric ere virtually non-existant 
[3.3 tee | Sven the molecular theorists have contented themselves 
with predictions of excess functions in the liquid phase at neg-. 
ligible pressures or when venturing into the field have confined 
themselves to predictions of relations for the critical constants, 
In the gas phase preoccupation has been With determining ane camae 
muon of state Ea and almost no experimental work has been done 
to determine excess properties of gases, Only very recently have 


badly needed measurements on dense gases been attempted [465,47]. 
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It had been generally felt | | AOE OB lels molecular theory 
would stand up well in the estimation of excess functions of mix- 
mares of condensed gases, that is,sliquid mixtures of simolesnol 
ecules at temperatures near the boiling points of the mixtures, 
moechne past decade, however, experimental work by Beenakker ecusal 
[2,3,4| and Staveley et al 6,7 | have shown serious inconsisten- 
mies 2nd plaguing weaknesses in the molecular methods. They have 
‘found that values of excess energy calculated by tne notecuia 
theory of Prigogine et al | 4-3 | are strongly dependent on both 
Bae set of potential parameters chosen and on the pure component 
femeceved 2s a "reference" substance, 

The most probable cause of feaeee OF the eo ie "Cons ir= 
dered to be the Lennard-Jones potential | 6,7,48,49, 50 | particular- 
ly at the close molecular distances involved in liguids. Still 
Peouner c&@use may lie in the een of the interactions rar 
emeter combining rules, At any rate, the spectacular differences 
in the value of excess functions calculated from molecular stiecons, 
from those measured experimentally is enough to precludeat this time 


Sonstacration of values formulated by .this means, 
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Salorimetry 


The measurement of the excess heat is of Such importance as 
to deserve a separate section reviewing devices and methods pres- 
ently available for its deternination, The two principal instru- 
mame are the adiabatic and isothermal calorimeters,  Unesacdi2= 
Baoic calorimeter reflects positive. excess heat, if tne tempera. 
ture S2OpsS On Maxine two pure components im 4@n adiabat iemee.. sae 
fone if On mixing or afterwards heat is added to cancel out the 
termperature drop, the amount of excess heat can be measured, 
Heat added during mixing reduces corrections necessary for heat 
transfer to the environment. A Similar determination may be made 
fOr negative excess near Witn its indicative temperavice sere 

The Sdtieiatie calorimeter has been poltnercdy py a Variety sone 
problems including imperfect heat isolation, incomplete mixing, 
inaccurate measurement devices, and svaceS in which vapor can 
form, The last problem is complicated by the need to allow for 
volume changes on mixing. . 

Peet nermal calorimeters have the advantage Of “noOG. Vea rae. 
Merrecctions for heat losses or gains and are able to measure a 
Meeaper Cemperature range BeRuRATeIGr: 

era Such devices are mentioned in the literature and Some 
of the more important are listed here, some early calorimetry 
feo Cemperatures of about 100°K Was fone by Rowlinson & Malcolm 
[2]. A versatile instrument which can be used either isothermally 
or adiabatically has been developed by ictnel Ser BE A glass, 
highly accurate, apparatus is described by Larkin & McGlashan [ 3]. 


mich has eliminated vapor spaces, provided conplete mixing and 





ellowed for exvansion and contraction throug] cenmpencaviaG 
mstcon, Likewise, FMrazek &« Van Ness 2: | Tie aegis equipment con- 
meol tne volume with close clearance pistons, Avot aSeaemel con 
trol on volume change is the metal aneroid diaphragms of the de- 
meee Of Brown & Fock 5]. Work has been done with Bunsen iso- 
thermal calorimeters by some investigators 16,7]. Two recent 


calorimeters were specifically designed for measurements of con- 


densed gas mixtures [8,9]. 
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Latent Heat 


Meemiacent heats of vaporization for mixtures are covered by a 
meititude of definitions [1,2]. The ost obvious are the mMeatcs 
meerired fOr vavorization or condensation at constant pressure, 
memevant tCennerature or constant composition. 

One method of calculating the latent heat is by employing 


wareacvions of the Clausius-Clapeyron relation [1,3-8]. One of 


LA 





these variations expresses the differential heat of vaporization 
Ya) yay ae | V il Zl We iN 
fow,,6P - (A\H,/T)aT - (x - x )(0°G*/Qx “Jp pdx = 0 (50) 


for a small amount of vavor phase out of or into a large amount 


meer guid, Similarly, the differential heat of condensacwene. 


Av, ar - (An, )/t)at ~ (xt - xY)(O2GY/ x72), pax” me (51) 


for a small amount of liquid phase out of or into a large amount 


of vapor, The volume change is given by 


Ay, =v’ -v> ~ (x" - x) OV OX Vp a (52) 


and the heats are 
V 1 V ai mi el 
BH, = ee Toes = CO Hn / Ox a (53) 


and represent the heat required at constant temperature and pres- 
sure to obtain one mole of vapdor phase from a large quantity of 
HyqQuid or to condense one mole of vapor into a large aceune on 
meed, Similar exvressions apply for A\V,, and {\H,_. Even with 
this seeninsly straigntforward Appooncls there is no Simple method 
Meee ion except with assumptions whicn ignore the effects of 
non ideality (see Graphical Methods). 

Corresponding States theory yields an imposing array of ests 
mation techniques for pure component latent heats 7]. Calcula- 
tlons were made by one of these methods lo], alleged to be among 
the most accurate and simple to apply, by adapting it to mixtures 
by use of the pseudocritical constants. This empirically derived 


relation can be expressed 





(Ata), = [Tp(5 108 Pon - 2.27)] /(0.93 - T,) (544) 


where Ty is the normal boiling point and AVN, is the isotherral 
earferential heat of vaporization. The results, caleutatedvat 

one atm, were less than satisfactory as can be seen in Fig, 31. 
ieract, because of the large latent heats 4S compared to tne ex— 
cess heats involved, it is a more accurate estimate simply to use 
the mole fraction avefage of the pure component latent heats, 
peor culariy at low pressurés, 

One experimental determination of the latent heats of No-O9 
mixtures exist 20}, a notable number when so few for any binary 
mixture can be found, Danats measurements at atmospheric pressure 
are plotted in Fig. 31. They represent the integral heat of vap- 
Meeation at constant pressure, that is, the heat required to 
@nenge one mole of liquid to the vapor state under conditions of 
constant pressure and composition, 

Me matter now tne heat of vaporization is determined, Nnow- 


Peer, it is only as useful as the isotherm or phase boundary from 


Which it will be based, 
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Graphical Methods 


(fer cem Test graprrcam tectrniques for determining the thermo= 


dynamic properties of N5~0 mixtures were devised at the Kamer- 


2 
lingh Onnes Laboratory at the University of Leiden Fra? Keesom 
& Tuyn developed a general construction for H,S,x Surfaces of ihe 
mixture in the two-phase region at atmospheric pressure, The 
property of Pics @maxture that is the key to tne geometric method 
is the transverse fold in the H,S,x-surface representing the wet 
vapor region. A tangent plane across the fold contacts the points 
of equilibrium or coexistant phases. These points were called 
Memmodes; the isotnerm joining them, the nodal line; and the 
Moeus of these points, the connodal curves ianouens liquid 
curves). | 
Basically, their procedure was to locate known points, such 
memeune bDOoliling point of oxygen on an H-S diagram, and with the 
aid of pure component heats of vaporization at l' atm, the defi- 
metion of enthalpy Mime myropy In terms of specific neat, ideal 
Mixture rules, and Dodge and Dunbar equilibrium data | 3]. Dio 
the vapor line which is then easily projected into the S-x plane, 
Using the previously mentioned princivle of tansent planes. 


the liquid line was determined in the H-S and S-x planes by 





Prapiically ee the line representing the plane tangent 
mo the connodal points, assisted by Dana's heats of vaporization 
at constant composition Gr The projection in the H-x plane was 
then easily established. 

The procedure then recommends extrapolation in the H-S plane 
Beechne Lines of constant composition to higher pressure using 


Second Law relationsnips and the virial equation of state, ie, 


dH = V = lV Bir ap and 
ox 
(550) 
fee= -(OV/ OT)p oP 
which become with the application of the equation of state 
Pt’, - P,)(B - TdB/dT) and 
(56) 
S=R ln P,/P, - (Py - P,) dB/aT 


The ideal mixing rules locate the isocomposition lines between 
those of the pure components, and the equilibrium data.fixes the 
end points at the vapor line. The form of the liquid line is 
then taken to be the same as that for one atm. 

Lerberghe 5 | proposed a variation on this method of his 
colleagues which requires only the use of the T-x equilibrium 
diagrams, From thermodynamics and analystical geometry, expres- 
Sions were develoved for the enthalpy and entropy difference be- 


tween tne connodal points. 


(N 4+ MT) (xY - xt) 


(57) 
ly 3 


ae 

< 
1 
26; 
i 


and 


(L + M)(x” - x 


7) 

< 
! 
op) 
It 





LO 


V V 
where L = SH aa S3 + R in [(1 = XK wy 
ie -RT(1/(1-x") + 1/(x")] ax" /ar (58) 
Vy V 
N = Ha a He 


Where the A and B are the pure component properties at the pres- 
Sure and temperature of the mixture, extrapolated as necessary | 
into the wet vapor re@ilon- 

Thus by fixing the vapor line as was done by Keesom and Tuyn, 
Lerberghe'ts method locates the liquid-line with the help of the 
equilibrium data, Note that the requirement exists to perform 
me difficult graphical tangent determination for ax /aT Or EEO 
@erive an equation for the vapor line on the T-x diagram by curve 
fitting and then taking its derivative analytically. Lerberghe's 
graphical measurements are reflected in Fig. 31 for the liquid 
mre At one atm. 

Several years later Bosnjakovic [6 | nicely summarized vari- 
ous grapnical procedures. His metnod She aoeee that in most cases 
the excess functions in the vapor phase can be neglected, at 
least in moderate pressures, and the isotherms in the H-x and S-x 
diagrams may be drawn according to the ideal mixture rules. Then, 
if tne temperature-vapor composition relationsnip is known from 
Pomelibrium data, the vapor line is also determined in both 
planes, 

In the liquid region, if the excess properties are known, 
then the isotherms may be drawn and the liquid line again deter- 
mined as above for the vapor line. If the excess heat 2 aoe, 


known, some estimation of the latent heat of vaporization must 





be made. The Clausius-Clapveyron relation (Sq.-51) makes a good 
meartcing, point. From the definition of Gibbs rece ecner. wr ie 


Becona derivative is given by 
2 2 2. Z Z 2 
(CO C/OX Vp m= CO H/OX Vp ne TOO S/O X) 5 a (59) 
Then for an ideal vapor mixture 
(O°G"/Qx"2),, = RE/x’ (1-x”) (60) 


Measurement of the slope of the T-x diagram graphically or analy- 
Gecally yields Opa Wa. Enough information is then available 
to solve Eq. (51) for eines at constant pressure and temperature, 
memevapor isotherm is then extended to intersect x* and (\ H, 

ms plotted from that point (Fig. 5). Repeating this process 

Pith msounerms Of the wet vapor region yields the liquid line 


(Fig. 31). Obviously the method is sensitive to the behavior of 


the two differential terms of Eq. (51). 


 . 


EE Auiv/? 


vapor line 





isotherm 





A ye B 


Figure 5 Figure 6 





In the case of the S-x diagram, after tne isotherms have been 
drawn according to tne ideal mixing rules and the vapor line de- 
termined from equilibrium data, the liquid line may be determined 


from the relation 


Uh 


a, /T ae Sen er Ome x )(OS /Qx), q (61) 


Tne value of AEE A can be obtained from the H-x diagram and 
mmetcea from the intersection of the tangent to the vapor linevar 
x” with the Sqr pr iumetiquid composut ten, = (Pic, 260s 

A secondary method of constructing the S-x diagram when the 
H-x diagram is available utilizes an integration procedure, which, 
mMeediess to say, is less prone to graphical error than the dif- 
ferentiation technique. This method uses the relation that, for 
a reversible vaporization process at constant pressure and con- 


Stamb liquid comvosition, 


Vv 


TY 
tt : 
.m - S Jean = fan (aH) 1 (62) 


ace 
mMne right hand side of the equation may be easily integrated 


graphically using any of the several numerical integration rules 


and the difference measured from the known value of eu 
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The Thermodynamic Properties 


The thermodynamic properties of binary mixtures of nitrogen 
and oxygen considered to be of the most value for use in engineer- 
ing calculations are those of enthalpy and entropy. A useful 
format for the presentation of these properties is the graphical 
eeeoplay, in particular the enthalpy—-composition and entropy— 
composition diagrams at constant pressure. Accordingly Figs. 
7-30 nave been prepared for mixtures at pressures of 1, 5, 10, 
and 20 atm. Enthalpy and entropy have been represented as func- 
meons Of composition in the liquid line region, the vapor line 
meeron, and the wet vapor region. 

Values for the pure component oxygen were taxen from Stew-— 
art's machine tabulation [2 | based on a complex empirical expan- 
mlomee tnose for nitrogen came from the work of share laine [2 | 
whose machine tabulation is likewise based on @n empirical equa-~ 
mvom of state, | 

The molecular weight of nitrogen was taken as 28,016 
grams/mole [2 | and one calorie as 4.184 abs. joules, by which 
conversion factors the tabulated values of Strobridge were con- 
verted to abs. joules/mole. 

for tne sake of graphical convenience the propertics of 
Peeerosen and oxysmen were equalized at 91°% at l atm, creating an 
oxygen enthalpy scale factor of 4245.3 j/mole and a nitrogen en- 


tropy scale factor of 81.12 j/mole. The values of entropy are 





dy 


relative only, Since the nitrogen scale is nov cased sone oe 
tropy at absolute zero as is that of oxygen. Urroperty oe ee 
the pure components within the two phase region have been linear- 
ly Bera noaaced: 

Values of the entropy of mixing (Gibbs paradox) are derived 
meom tne relation 


ee R _ nx, (63) 
A. 


where R is the universal gas constant equal to 8.20574 x aioe 


L-atm/mole-°K 1]. 

Isotherms in the superheated vapor region were derived as 
mole fraction averages of the values for the pure components, ie, 
meemecne ideal mixture enlee, The excess heat and entropy func- 
tions for the liquid isotherms were derived by the author from 
the data of several investigators (see Discussion of Results). 
fuge Vapor and liquid lines were plotted using the isotherms and 


the temperature-composition data of Dodge and Dunbar [3]. 


References for Thermodynamic Proverties 
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DISCUSSION OF aaa SotiibEs 

The cnoice of the metnod of determination and the form of 
presentation of the results was arrived at after consideration of 
all the previously discussed techniques. No attempt was made to 
formulate a virial expansion an other empirical equation of 
state, because of a general lack of the PVT data required for 
such a formulation, A corresponding states theory vesrimeriem 
technique, uSing pseudocritical constants and the enthalpy devi-e 
ation tables of Lydersen, was conducted, with the results already 
mentioned, The negative excess heats calculated by this method 
are thermodynamically inconsistent, and this procedure was re- 
jected. Using pseudocritical constants, a modified relation for 
the latent heat was calculated, which when plotted produced the 
least agreeable of all liquid lines determined (Fig. 31). The 
purely graphical approach of Lerberghe vroduced better results, 
Seeeoy AL hign molar concentrations, which plotted extremely 
erratically.. Both a mole fraction average of the heats of vap- 
@riZation of the pure components and a mole fraction average of 
Baewpure component boiling points produced better results in 
that order, than did the comvlex estimation techniques mentioned 
above, A plot of the results obtained by these methods at 1 atm 
‘together with the method finally chosen is given in Fig. 3l. 

ia the end ‘it be decided that the calculatloneeiecess 
minmeelOns provided the best determination of the liquid line, and 
therefore of the isotherms in the liquid region. The liquid 


line determined by this method shows very good correlation with 


the experimental evidence of Dana, as can be seen in Fig. 31. 
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Figure 31. Liquid line for N5-O5 mixtures at atmospheric pressure. 


A presentation of the excess Gibbs free energy as calcula- 
ted by previous investigators is shown in Fig. 32. The excess 


Gibbs free energy is related to the other excess functions by 


meec 4s’ - Gc . 7 O6°/D7), (46) 


Mmemtoecliear from Fie. 32 that the slope of ea with respect to 


termperature is very linear in the low pressure region and has 
been measured as 0.08 cal/mole °K. ‘ith this slope the excess 


heat of the liquid was calculated to be 16.0 cal/mole at equimolar 





fF I 


composition and low pressure. Only one experimental calculation 
er excess heat for N5-O5 liquid mixtures has ever been made [1]. 
These measurements by Knobler et ai of the excess heat at nk 
Pee OLOtted as Fig, 33. These results can be expressed Dye vic 


fitted polynomial approximation [2]. 


a. =. 
H’ = XNo Xoo een 10(2x, 


ae 1) + BI ae = 1)? |ca1/mole (64 ) 


which yields an excess heat of 10.5 cal/mole at x=0.5. A concur- 


nw 
rent calculation of G was made and can be given by 


E : 
G = xno XQ5 | 40 A 4(2xXo5 ~ 1) + 4 (2x0, - 1) |ca1/mole. (65) 


This gives an excess Gibbs free energy of 10.0 cal/mole at x=0.5 
or a value nearly the same as that of the excess enthalpy. Since 
the temperature derivative of ae is negative, then Hee. and 
from the weight of evidence by considerably more than 0.5 cal/mole. 
‘It was decided that 16.0 cal/mole was a reasonable value for the 
excess heat at equimolar composition in keeping with Fig. 32, 

and the sood fit obtained in plotting the liquid line? Knebler = 
data was then scaled up to this result. It should be noted that 
in agreement with prediction [3] SINS GOSS heat is asymmetric 
With respect to the composition, while the curve of excess free 
pwemey 1S Symmetrical. Values of the excess entropy were then 


calculated from 2a. (46). 

When the calculated values of excess heat were applied to 
isotherms in the liauid region at one atmosphere, the resulting 
liquid line showed excellent agreement with the experimental 


Gace of Dana, It was decided, therefore, that this technique 
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Figure 32. Excess Gibbs free energy as a function of temperature 
for equimolar liquid mixtures of No-Oo. See Appendix C for data 


sources. 
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Figure 33. Excess heat and excess Gibbs free energy as a function 
of composition from the data of Knobler et al. HE fitted to re- 


vised equimolar avlue of this thesis. 





Ge 


at atmospheric pressure was the best available, 

The integration method of Bosnjakovie was used to calculate 
es ~ a ee The liquid line thus determined for the S<xsdace 
gram was in excellent agreement with that calculated from values 
of the excess entropy. 

In the vapor phase it had been decided that most authors 
generally agreed that there was negligible heat involved with 
mixing at low pressures including atmospheric, so that the best 
approximation of isotherms was a Simple mole fraction average 
corresponding to the mixing rules for an ideal mixture ees 
Cludes the entropy of mixing for aa diagram vapor line, 

The major problem then, concerned pressures greater than 
atmospheric, that strangely ignored and abandoned region where 
investigators have feared to tread, It would be foolhardy and 
certainly not justifiable to make unwarranted assumptions about 
conditions at higher temperatures where pressures exceed: several 
atmospheres, where the third virial coefficient arc comoressi- 
pility of the liquid, and the volume change on mixing are all 
important factors, 

For the best available presentation it was decided that the 
excess functions of heat and entropy already calculated would be 
used to determine the isotherms of the liquid and mole fraction 
averaging the vapor isotherms at pressures and temperatures be- 
Pow the critical. This decision is not wholly without justifica- 
tion. It is not unreasonable to assune that the excess thermo- 
dynamic properties of the liquid are relatively independent of 
pressure fel, vrincipally because of the relatively small volume 


changes in the liquid state, Additionally, Fig. 32 gives some 





a 


small nint that deviations at the higne* temperatures associated 
with the higher vapor pressures @re not so great from the more 
mecurately calculated linear results Me lower temperatures, 

Tne incompressibility reasoning does not apply so well in 
The dense gas region, but here too it is felt that the ideality 
Se mixture assumption is valid enough when compared to other 
available data, The values of several investigators for air 
[5.6.7.8 | were compared at the pressures in question at several 
{sotherms near the dew point and compared to the 79%mole No mix- 
ture of nitrogen and oxygen isotherms calculated by mole fraction 
averaging. Although the deviations varied widely, a positive 
excess heat with increasing pressure might seem to be indicated, 
the effects of argon in the mixture and the particular investi- 
gations notwitnstanding. The deviations are not large enougn nor 
their pattern clear enougn to justify guesswork at possible val- 


ues of excess heat, 
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CONCLUSIONS 

ice avemuron. Of 9 unis thesis was to survey the literature 
with regard to the thermodynamic properties of binary mixtures of 
mMecrogen and oxygen in the two-pnase region and to analyze the 
resulting data. From,the sizeable collection of available print- 
ed matter which concerns itself directly and indirectly with the 
Subject, it would seem an easy matter to make an analytical de- 
termination of the common thermodynamic properties of such rela- 
tively simple mixtures, A wide variety of estimation techniques 
do, in fact, exist, and while those have been presented here, 
most, if’ not all, of them are either unsatisfactory or at least 
only partially useful. 

Specifically, the virial expansion and other equations of 
meawe With their associated derivatives would be most satisfac- 
ory if experimental PVT data existed with which to determine 
moar coefficients. Molecular CNECOry, Simidgariy, nas s<reatl poe 
Seed particularly for liquids and-dense gases, but suffers 
meom tne uncertainty involved in its, basic parameters, The tech- 
menagues involving pseudocritical constants give good rough esti- 
mates, but do not vrovide the accuracy required for even engi-e 
Meerans calculations. The excess functions provide a rapid and 
relatively easy method for determining the thermodynamic proper- 
ties. Presently available data on these give good results at 
pressures near atmospheric, | 

All of the aforementioned methods have one thing in common. 


None of them have been developed for use in the high pressure 





area of the two-phase region. Where no one has gone before, how- 
ever, any determination, no matter how simplified, idealizeqdyec: 
approximated it may be, 1s better than none at all. Therefore it 
Ba decided that by taking the excess functions independent of 
pressure in the liquid region and using ideal mixture rules in one 
vapor phase, satisfactory thermodynamic diagrams could be con-~ 


structed, 


RECOMMENDATIONS 

Under this heading of Recommendations Should be listed the 
solutions to the problems posed by the Results. At the risk of 
redundancy and speaking only for the mixtures in question these 
inadequacies are repeated here, 

A Single equation of state for all mixtures would be the 
most satisfactory solution as it would provide in one analytical 
expression and its derivatives the measureable and derived prop-~ 
meewtes fOr the infinite number of mixtures. A .solution of this 
type would require an extensive experimental Beorren to deter- 
mine a sufficient number of PVT data to permit the evaluation 
of appropriate coefficients. 

Another satisfactory program would be to make direct measure-~ 
ments of the excess heat faneeien using one of the calorimeters 
mentioned previously. This is the only satisfactory means of 
determining this property and at high presSure the only means. 
Then since the vapor pressure method of determining G™ is only 
satisfactory at essentially senoserente pressure, the excess en-~ 


tropy function could be determined by the integration procedure 





discussed earlier under Graphical Methods, A program of this type 
would not be all inclusive in its results and would necessitate 
taking measurements at sufficiently close intervals of temperature 
to permit reasonable interpolation of tabulated data or accurate 


drafting of thermodynamic diagrams, 
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Appendix A 


1. Temperature-composition 
2. Temperature-composition 
3. Temperature-composition 


4, MTemperature-composition 


elalieVescreil. Ml rela - 
diagram, 5 atm, 
diagram, 10 atn, 


diagram, 20 atm, 
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APPENDIX B 
Barker's method: 
The most widely noted method of determining ee from measured iso- 
thermal values of total vapor pressure and liquid composition is 
wnat set out by Barker, 
The excess chemical potential is given by Eq. (47) which can 


be written 


0: ae a V i Ih 
My. =e 1 No = RT In (Pxy5/Py OX) + (Vivo - Byo)(Pyo - P) 
GaN 
+ P On509(%O5) (66) 
With a similar equation for | oo 
It follows that 
PS VP uo : F 0% 02 (67) 


where Pyo and Fes represent partial pressures if the liquid is an 
‘ideal solution but the vapor is a non ideal mixture. They are 


Biven by 


= It 1 : vay 2 
Py = XNoPyp CXP (Vip - Byo)(P - Pyn/RT - PO ng09(XOq) “/RT (68) 


with a similar expression for P.. and where 


02 
ONn202 a 2B 505 = BN» - Boo (69) 
lea 
me 1 pe Life . =e Le | Lr gucelemay 2 ences 
= Xo no + XooMOos = *No Xoo a + b(xy,, X55) + e(Xn.-XG5) 


(70) 


then 





in No = A ENS + B MN» + C "No ; (7a) 


and similarly for lin Too 


The values of 1, m and n are expressed as 


1 aes ee 
ae "05) ’ NO os) 
le? il 62 1 
ih. 2 1 ee ee 2 1 1 eZ 
Ny o=(Xo,) 1 - 8xNo + AGE No n=(%y,) 1 - 8x0, + LM) 


and A= a/RT, B=b/RT, C=c/RT, 


As a first approximation the solution is assumed to be regular, 
Baeimecnat B=C=0O, By neglecting corrections for the vapor phase 
non ideality, A is set equal to 41n(2P VA + Py es 


ie 
Es and Po are then calculated for the various experimental 
2 2 
liquid compositions from Eqs. (68), These are used to calculate 


Tos? Ux, and Poi, from Eqs. (71) & (67) which in turn yield 


the pressure residual 


R = ot = EEA (73) 


and the derivatives 


= Pp P ly 
aP/dA ba b wo? + lo, F 00° 0. (74) 


with similar expressions for dP/dB and dP/dC using m and n re- 
spectively. 

Changes in A, B & C to most nearly reduce the pressure residuals 
to Zero are determined by fitting by least squares to 


(4P/aa)OA + (AP/a3)OB + (aP/aC) Oc = R (75) 
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or by solving simultaneously equations (ot st lesnerm 


Sa) (ap/aa)? + Sp) (ap/aa)(ap/aa) i Sc) (aP/aA)(aP/dC) = 
) R(aP/aA) (76) 


with similar expressions for ) 2 (aP/aB) and ) R(aP/aC). 











Appendix C 
Selected Data 


ie 2Quilibrium Data 


P T xl xv Ht HV sl sv 
Atm °K mole No Zmole N5 j/mole j/mole j/mole °K j/mole °K 
i Bo BE IONS 100 823.7 6793.5 92.84 165.00 
78 89.8 97.3 722.0.6790.00 95.62 166.23 
80 64.0 BA W722. Ore ae 98.00 Date) 5 ils 
82 Hh , 3 ee Ase PEM olsiclowne SGm27 LAs 
Buy 29.6 63.8 160.0 6607@ O7 58 r7 3 Oe 
86 1880 47.6 70.0 G5vPmO 97808 17ECO 
88 8.3 2; (0) 8.0 6486.0 95.95 173.88 
90 un ae =P ONO 964 Come 94.50 170370 
90.18 0 0 23270 6 oom 93.94 169%.63 

5 94.18 100 100 Pele 5 6672 43m eal 155.84 
96 79.3 91.0 16301076755. Ome Oc 2e 159.05 

98 Cll 80.4 149220 6GS00nM0un Loo .29 ieedes 6) 5) 

100 4o.? i 1337.0 6865 70m 109,60 163.09 
iO? poeo Ce a7 1260 6O2RmOmeslOo. 26 164,32 
104 220 2 1171.0 6994.0 108.58 164.80 
106 2.3 2Oe6 1108.0 7053.0 107.45 164.58 
108 B56 oa | 1O57 20s 7 Memon. 105.76 162.90 
108.99 ©) ) WS), ash ae THO EG, 160.63 

meme O3.95 100 100 20.2 670205081 10, 34 151,36 
104 99.3 99.5 2433,0 6702508 210,57 lieaiag SiS 
106 Von 89.0 226550) 67 Soom 27 .155.12 
108 63.6 eee 2137.0 6855.0 115.38 157.24 
10 49.6 68.0 2026.0 CO2ZC ROIS oS 158.88 
BZ SHAS one 1939.0 6994.0 115.49 159295 
mie 26,3 43,4 1858.0 7060.0 114,89 160.55 
116 Ie 2 ZO AL L7G. 71 Zee neld 34079 160.50 
118 OAS ISAs 172,007) Sono emer eae 159.30 
19.86 0 0 1681.7 7245.1 110.09 156.50 
Poeeli5.82 100 HOO 3321.9 65152 kl? .96 145.53 
116 Bie 2 99.0 3315.0 6529.0 118.48 146,00 

ma 8 82.3 S18) 5 IL S77. 0 Gol Ome meet, 52 149.78 
120 68.8 79.4 3072.0 67 5ememel2c.79 151.98 
ieee 56.1 69,2 2072.0 6S jomomel2 3. 09 153.70 
124 ee 58.0 2373205 CO NeGmomere2 9S 154.91 
126 33.0 46.0 2725. 0 OOO mOmE 22 . 6 155.62 
128 228 Bon7 2767.0 7O05emOmmeIL 2). 53 WSS 577 
130 Insial 20n7 PEG. 0 72cm 20.13 155.23 
132 bed fares. 2543.0 7190.0 118.13 LSS 
133.03 ©) 0 254250 20 woer b16.ol iS eo 





2. Fure Component Data 


Pp di HO Hy S Sx 
Atm OK sy One {i/o j/mols ox symone OK 
1 65 =) 347 108.4 76.95 8277 

70 ~1115.3 397.6 80.38 87.05 

75 =O Glee 687.0 83,88 91.05 

80 -594 ,3 6486, 3 7 soe INES, CF 

82 ~485,3 6546,8 Senor WS 777 

84 -375.9 6607.3 89.99 167 aoe 

86 =265e2 6667.8 91.30 lesneu 

88 = Sh 6728.0 92.57 168.90 

90 ~42,3 678eme 93,83 169.56 

95 6940.0 6938.4 7 ieee 72 © 

100 7001. 7088.1 MABE U7 

5 80 SSA, 5 983.1 87.28 94,68 
. 85 ~313.5 1270.8 90.60 98.18 
90 =i 356 Soll © OBO 101.52 

96 300.9 673700 97.40 BY) Gul 

100 526.4 6876.3 99.70 Sey oH 

104 752.9 7012, eo 159.26 

108 981.0 7145.8 Oe 7 160.52 

Ae 7179.5 HP 160,94 Ges 

115 7349.4 7374.9 162,45 162.56 

IAG, (ones, 5) MSs) 0 0 WES 57 163.93 

10 100 Sone Pe a6 99.63 107.79 
104 759.8 6704.5 101,84 151.38 

106 C7362 6790.2 OZ me 152.20 

108 937.1 6872.9 103.99 ‘152.95 

110 1101.6 6953.3 105.04 153471 

We AMS 9, (AOL 7 106,07 154,41 

114 1333.1 MOS. S ILO Le) LES 50S 

LS 1450.5 TVG 326 Oee lez 155.75 

118 1569.5 72 5a 109.14 156. 37 

20 C2507 (a DUSY SiO 5 156.99 

20 105 830.3 2508.0 e223 110.59 
LUG jbl) 5 2856.0 104,86 113,84 

115 1400.8 3249.0 Ogee II Abe 

118 1576.5 6663.6 108.92 146.89 

1.22 1816.7 6894.7 EO. 2 148,72 

126 2066,4 7096.5 112.94 150.35 

130 2329 a7 7280.5 ime Pnele 151,80 

135 7olaine ILOL .8 152.46 153.40 

140 (iow (255: 7697.1 154,21 154,89 

145 07 ee 7890.4 ie Se7 Gece 





3. Deviations 


D6 


mole N 
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4, Summary of equimolar en 


allen = 3\S) 
notation 


Oy me 
63.14 10.88 
65 Hoes 
70 10,28 
Gs) 9.74 
75.16 9.65 
7765 10.75 
ae 10.0 
Bono? 9,2 
90,0 =« 9.6 
91.0 (Ss 
80 845 
90 Soe 

100 9.75 
110 9.08 
120 7.64 


Source notation refers to "References for Excess Functions" 
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